Tunable High-Frequency Gravitational- Wave Detection with optically-levitated 

sensors 
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We propose a tunable resonant sensor to detect gravitational waves in the frequency range of 30 - 
300 kHz using optically trapped and cooled dielectric microspheres or micro-discs. The technique we 
describe can exceed the sensitivity of laser-based gravitational wave observatories in this frequency 
range by 1 - 3 orders of magnitude, using an instrument of only a few percent of their size. Such 
a device extends the search volume for 100 kHz gravitational wave sources by more than 10 6 , and 
could detect monochromatic gravitational radiation from the annihilation of QCD axions in the 
cloud they form around stellar mass black holes within our galaxy due to the superradiance effect. 



PACS numbers: 04.80. Nn,95.55.Ym,14.80.Va 

Introduction. Over the past 40 years optical trapping 
of dielectric objects, both macroscopic and atomic, has 
made a profound impact in a wide range of fields rang- 
ing from fundamental physics to the life sciences. First 
studied by Ashkin and coworkers optically trapped 
dielectrics in ultra-high vacuum become well decoupled 
from their room temperature environment 043- Recent 
work suggests that the center of mass motion of such 
levitated objects can attain mechanical quality factors in 
excess of 10 12 , while internal vibrational modes are com- 
pletely decoupled. This remarkable decoupling can be 
harnessed for cooling the center of mass motion of such 
objects to the quantum ground state @-@|. These sys- 
tems also have been considered in the context of reach- 
ing and exceeding the standard quantum limit of posi- 
tion measurement [9J. In addition, these techniques en- 
able ultra-high sensitivity force detection [Iol - fl2| and ex- 
tend the study of quantum coherence to the mesoscopic 
regime. 

In this paper, we study how nano- and micro-scale sen- 
sors trapped inside a medium-finesse optical cavity can be 
used to detect high frequency gravitational wave (GW) 
radiation. While there has been convincing indirect ev- 
idence for the existence of GWs [13l | . their direct detec- 
tion has yet to be demonstrated. Such a detection is 
highly likely in the next decade with the new generation 
of laser-interferometer gravitational wave observatories 
and will launch the field of gravitational wave 
astronomy. While these detectors have been optimized 
in the frequency band of 10-10000 Hz, their sensitivity 
decreases at higher frequency due to photon shot noise. 

We propose an alternative form of detector for im- 
proved sensitivity in the frequency range of 30-300 kHz 
that does not rely on a shot-noise limited displacement 
measurement of test mass mirrors, but rather depends 
on a precision force measurement on the resonant har- 
monically trapped sensor. The method we describe can 



[vwwv] 



FIG. 1: A dielectic nanosphere or microdisc is optically 
trapped in an anti-node of a cavity of length £ m at position 
x s . A second light field (not shown) is used to cool and read 
out the position of the levitated object. A passing gravita- 
tional wave at frequency u) gw displaces the sensor from its 
equilibrium position in the optical trap and imparts a force 
as described in text. The resulting displacement is resonantly 
enhanced when u) gw coincides with the trap frequency loq. 



yield sensitivities improved by 1-3 orders of magnitude in 
this frequency band, and thus extend the effective search 
volume for sources at around 100 kHz by over 10 6 . 

Finally, we discuss GW sources at high frequencies. We 
focus on GW signals from the effects of the QCD axion on 
stellar mass black holes (BHs) through BH superradiance 
[l9l | . This novel signal comes from axion annihilation to 
gravitons, is monochromatic, long-lived and extremely 
different from all known astrophysical sources. 

Experimental Setup. The apparatus is schematically 
shown in Fig. [TJ We consider a dielectric sphere or mi- 
crodisc with dielectric constant e, mass m, and density 
p optically trapped and cooled in a cavity of length £ m 
using two light fields of wavevector k t = 27r/At ra p and 
k c = 27r/A coo i, respectively. The sphere or disc, which 
we generally refer to as the sensor in the following dis- 
cussion, is levitated in an anti-node of the trapping light 
located near one of the end mirrors (e.g. within £ m /100). 

In the locally Lorentz (LL) frame with origin at the 
input mirror, a gravitational wave, h = h(t — f ), that 
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propagates in a direction perpendicular to the cavity axis 
will change the proper distance between the mirrors as 
well as the distance between the mirror and the dielectric 
sphere or disc (see Fig. [T]), respectively: 



Quantity nanosphere 



microdisc 



L m = (1 + -jh)l m , 
X s = (l + -h)x s . 



(1) 
(2) 



With l m and x s we denote the proper distances between 
the mirrors, and between the mirror and the sensor in 
the absence of the GW. For the GW frequencies we are 
considering, corrections of order h-r? — can be safely ne- 
glected. The electric field component of the trapping 
light with frequency u) = k t c in the cavity is the sum of 
two counter-propagating waves given by: 
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In the absence of a gravitational wave, the sensor 
is located in an antinode of the trapping field at po- 
sition x s — x m i n . In contrast to the sphere position 
x Sl the change in x m i n in the presence of the GW is 
not determined by the proper distance change between 
two massive particles but by the interference of the two 
counter-propagating waves in the cavity. The condition 
for the antinode in the absence of the GW is given by 
kt(( m — x m - ul ) = (n + l/2)ir for an integer n. By requir- 
ing the same condition is satisfied under the influence 
of the GW, we obtain the new coordinate position of 
the trap minimum: k t (L m - X^ nin ) = k t (i m - x min ). 
The shift of the intensity maximum is thus given by 
<5^min = h(lm ~ x s )h + 0(h 2 ). Taking into account the 
effect of the gravitational wave on the sphere or disc it- 
self we find that the sensor is shifted with respect to its 
equilibrium position to first order in h by: 
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Eq. H] presents us with an intuitive picture of the GW 
effect on the apparatus. If the sensor is located exactly in 
the middle of the cavity there is no effect from the GW, 
as if the sphere or disc was connected to the mirrors by 
equal length springs while the distance between the mir- 
rors periodically changes. When the sensor is positioned 
away from the cavity center, the GW drives the sensor in 
an oscillatory motion away from its minimum, as if the 
sensor is now connected to the mirrors with springs of 
unequal length. 

Treating the GW as coming from a monochromatic 
source of frequency w gw and amplitude h a , there is a 
oscillatory driving force on the sphere: 
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TABLE I: Expression for trapping and cooling parameters for 
nanospheres and microdiscs. 



When the trapping frequency of the sphere matches the 
frequency of the GW, the sphere will be resonantly ex- 
cited, allowing for the detection of gravitational radia- 
tion. In this way the device operates like a resonant-bar 
detector, except with a very wide tunability [20j |. 

For concreteness, we consider a cavity of length £ m = 
100 m, finesse J- = 200, and cavity mode waist w = 75 
jum. We consider a silica (e = 2,p = 2.3 kg/m 3 ) 
nanosphere with radius a — 150 nm or an apodized mi- 
crodisc with radius a — 75 /xm and thickness t = 0.5 
fim, respectively. In Table HI expressions for the trapping 
frequency uq, trap depth U, gas-damping coefficient 7 S , 
and scattered photon-recoil heating rate "f sc are given in 
terms of the average intensity over the disc or sphere sur- 
face It, the background gas pressure P and mean speed V, 
the volume of the microdisc or sphere V, and the cavity 
mode volume V c = ttw 2 L/A, following Refs. [ 5|.|2 ll|. Here 
Tdisc is the micro-disc limited cavity finesse [2f| and the 
optomechanical coupling of the cooling mode is g. 

The cavity is driven with a trapping laser of wavelength 
Atrap = 1-55 ^m and power P t = 280 mW or P t = 385 
mW for a sphere and disc, respectively, corresponding to 
an axial trap frequency of u!q/2tt = 100 kHz. The cool- 
ing light has wavelength ~ 1.55 /xm, frequency uj c , input 
power P c = 50 mW or I /iW for the sphere or disc, respec- 
tively, and an optimized red detuning of 5 = —ujq. The 
cavity decay rate is k — ttvq /J 7 , and vq is the free-spectral 
range. Additional cavity loss due to photon scattering is 
negligible: less than I0~ 3 k for our parameters. The mi- 
crosphere or disc absorbs optical power from both the 
trapping and cooling light in the cavity, which results in 
an increased internal temperature Ti nt . Assuming negli- 
gible cooling due to gas collisions, the absorbed power is 
re-radiated as blackbody radiation. Tj nt is listed in Table 
im for fused silica with dielectric response e = ei + it2, 
with ej = 2 and £2 = 1.0 x 10~ 7 as in Ref. [25(, and 
£bb = 0.1 as in Ref. for an environmental tempera- 
ture T ext = 300 K. We assume T- lni and Tqm are not sig- 
nificantly coupled over the time scale of the experimental 
measurements at P gas = 10~ n Torr. Other experimental 
parameters are shown in Table HT1 

The Gaussian profile of the trapping beam near the 
mode waist provides transverse confinement, with an os- 
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dilation frequency of ~ 320 Hz. The earth's gravitational 
field displaces the equilibrium position of the sphere or 
disc by ~ 2.2 /jm. Cooling of the transverse motion can 
be done with active feedback to modulate the power of 
a transverse cooling laser using the signal from a trans- 
verse position measurement, for example generated by 
measuring scattered light incident on a quadrant pho- 
todiode. We assume a modest cooling factor of < 100 
in the transverse directions to counteract the effects of 
recoil heating and localize the sensor. 

For detecting the position of the sensor, light from the 
cooling laser can be frequency shifted to be on-resonance 
with the cavity to maximize sensitivity. The phase of 
the detection light reflected from the cavity is modu- 
lated by the sensors motion through the optomechani- 
cal coupling duj c /dz = k c g. Photon shot-noise limits 
the minimum detectable phase shift to Stp s» l/(2y/l) 
where I = Pd/(tiui c ) [22j. The corresponding photon 
shot-noise limited displacement sensitivity is y S z (uj) — 



for an impedance matched cav- 
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ity. In our case, a detection power Pd = .1 mW or 50 
mW for a disc or sphere, respectively, corresponds to 
yJS z {uo) = 4.1 x 10~ 16 m/VSz, and 1.8 x 10" 12 m/VSz. 
The thermally-driven resonant motion of the sensor is 
typically much greater (> 20 x) than the photon shot 
noise limited detection floor. We assume that substrate 
vibrational noise, electronics noise and laser noise can be 
controlled at a level comparable to the photon shot noise. 

The cooling serves to damp the Q factor to Q c g so 
that perturbations to the system ring down within rea- 
sonably short periods of time, to reduce the require- 
ment on the laser intensity stabilization, and to mitigate 
heating due to the recoil of trap laser photons. At the 
same time the mode temperature is reduced to T c g . The 
force sensitivity scales as ■uT e g/QeB- the minimum de- 
tectable force due to thermal noise at temperature T e g is 
-Pmin = \ Akk BT ctt b wriere j s the center-of-mass mode 
spring constant, and b is the bandwidth of the measure- 
ment. The thermal-noise limited minimum detectable 
strain due to a GW will be approximately 
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TABLE II: Experimental parameters for trapping and cooling 
a silica sphere with radius a = 150 nm in a 10 m cavity and 
for a microdisc in a 100 m cavity. 
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FIG. 2: (color online) Strain sensitivity for optically levitated 
micro-discs (black) or spheres (blue) for experimental param- 
eters described in the text. For comparision, also shown are 
the extrapolated LIGO and extrapolated predicted Advanced 
LIGO sensitivity in the frequency range of 10-300 kHz [2(|. 
The shaded region denotes predicted signals due to Black Hole 
superradiance as discussed in the text. 
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which describes the im- occurs. Here for a nanosphere the sensitivity scales as 



portance of photon recoil heating 7 SC and the efficiency 
of the cavity cooling. The factor R + , defined in Ref. 
Q can be minimized by going into the resolved side- 
band regime and can be generally neglected when com- 
pared with 7 SC . There are two general regimes of scaling, 
X « 1 and x >> 1. For x << 1, the effects of photon 
recoil do not significantly degrade the force sensitivity, 
and for a nanosphere, ft-ii m oc WQr 2 T~ 1 / 4 P~ 1 / 4 . In the 
regime \ » 1, photon recoil heating becomes signifi- 
cant, and damping without an equal amount of cooling 



hn m cx wo and is independent of r, T and P. The micro- 
disc geometry scatters much less light, as pointed out in 
Ref 2l| and recoil heating is significantly reduced. Also 
the micro-disc allows a larger mass to be trapped and lo- 
calized in a particular anti-node of the standing wave in 
the cavity. We assume a disc-limited cavity finesse of 10 5 , 
which is reasonable for an apodized disc [2l|. In Fig. [5] 
we plot the expected gravitational wave strain sensitivity 
using the selected experimental parameters. 

We expect a passive vibration isolation system similar 
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to that employed in LIGO would be sufficient for the pro- 
posed measurements, particularly since the effects of me- 
chanical vibration become less significant at the high op- 
crating frequencies. To properly distinguish gravitational 
waves from other disturbances, two perpendicular arms 
of an interferometer can be used, as in existing GW ob- 
servatories. Also, coincidence between several operating 
sensors can be used to discriminate other backgrounds. 
Motion in the end mirrors or their coatings due to ther- 
mal noise are not expected to significantly affect the dis- 
placement measurement at the high operating frequen- 
cies. Low frequency fluctuations or drift in the cavity 
length or laser frequency can modulate the intra-cavity 
intensity and therefore cause fluctuations Sluq in the trap- 
ping frequency. We require Sujq << wo/Qeff, which can 
be accomplished with active stabilization. 

For a sphere, additional motion (e.g. due to rotation) 
generally occurs at different frequencies and can be av- 
eraged out in a measurement. For a microdisc, although 
the restoring force of the laser tends to keep the disc 
upright, there can be torsional modes and axial rotation. 
The frequencies of these modes can be separated from the 
center of mass mode. Also, flexural modes within the disc 
can be neglected if their frequency is large compared to 
the variation of the effective trapping frequency due to 
laser intensity variation over the disc surface. We assume 
the disc can be treated as a rigid object for the proposed 
parameters, as the fundamental mode (2,0) of its defor- 
mation has a frequency of 126 kHz, which is larger than 
the variation in trapping frequency over the extent of the 
disc for u!q/(2it) — 10 5 Hz. This variation can be further 
reduced by increasing the beam waist to 1.5a, which still 
allows J^isc ~ 10 5 [2l|. This approximation becomes less 
accurate at higher trapping frequencies. 

GW sources. The most interesting source of GWs 
in the high frequency regime arises from the effect of 
the QCD axion on astrophysical stellar mass black holes 
(BH) through the Penrose superradiance process [Tij . 
The QCD axion is a pseudo-Goldstone boson that nat- 
urally solves the strong CP problem in strong interac- 
tions and explains the smallness of the neutron's elec- 
tric dipole moment 27]. Non-perturbative QCD effects 
generate a cosine potential that gives the axion a mass: 
H a = 6 x lQ- w eV/c 2 (10 16 GeV//„) , where f a is the ax- 
ion decay constant. 

The Compton wavelength of the QCD axion with 
fa ^ 10 16 GeV matches the size of stellar mass BHs 
and allows for the axion to bind with the BH "nucleus" 
forming a gravitational atom in the sky. The occupation 
number of the levels that satisfy the superradiance con- 
dition grows exponentially extracting energy and angular 
momentum from the BH. At the end of the process there 
is a Bose Einstein axion condensate cloud surrounding 
the BH. Axions in this cloud produce gravitational radi- 
ation through transitions between different atomic levels 
and through annihilations of axions to one graviton. The 



latter process becomes possible since the axion is its own 
antiparticle and the BH "nucleus" makes sure that en- 
ergy and momentum are conserved, exactly as a positron 
can annihilate with an atomic electron to a single photon. 
For annihilations, the frequency of the produced GWs is 
given by twice the mass of the axion: 

/=145tHs (iiii^Y) (7) 

which lies in the optimal sensitivity range of our setup 
when f a is between 10 16 and 2 x 10 17 GeV. The signal is 
coherent, monochromatic and thus completely different 
from all ordinary astrophysical sources. 

The GW amplitude coming from annihilations is: 

*~»-(?K^)fe) o) 

where a — GnMbhUci and £ is the orbital quantum 
number of the super radiating level [l9(. The ratio j 
is constrained by the superradiance condition and it can 
acquire a maximum value of ~ 0.5. Finally, e is the frac- 
tion of the BH mass the axion cloud carries and it can be 
as high as 10 -3 . In Fig. [2] we compare the annihilation 
signal for a source 10 kpc away to the experimental sen- 
sitivity assuming 10 6 sec of integration time. The slope 
of the curve is determined by the maximum BH mass 
that can superradiate for the given axion mass and from 
Eq. [5] we infer that it decreases linearly with increasing 
frequency, i.e. increasing axion mass. The curve satu- 
rates when this maximum BH mass equals the smallest 
possible BH that can be formed through astrophysical 
processes, « 2M . For a GUT scale f a axion, a sig- 
nal coming from within our galaxy could be detected. 
The number of BHs within our galaxy in estimated to 
be 10 7 — 10 9 and the axion annihilation signal can last 
for a few weeks or up to a year, so it is likely that the 
proposed setup has a good chance of detecting it. 

Astrophysical gravity wave sources have a natural up- 
per bound on the GW frequency they can produce. 
This is determined by the lower bound on the black 
hole mass produced through ordinary stellar dynamics: 
fmax ~ c /GjfMBHmini where Gjv is Newtons constant 
and MsHmin is the minimum black hole mass. This 
places fmax ~ 30 kHz for a one solar mass BH which 
is on the edge of our optimal experimental sensitivity. 

Discussion. The method we have described promises 
to be the most sensitive approach for detecting GWs in 
the frequency range over 10 kHz. Although few astro- 
physical sources are likely to exist at such high frequen- 
cies, there can be a variety of sources associated with the 
early Universe and Beyond the Standard Model. Such a 
source is the well-motivated QCD axion. 

The minimal requirements on vibration isolation at 
high frequencies relax limitations on the geographic lo- 
cation of the experiment. Several beads or discs can be 
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simultaneously trapped in different anti-nodes of the cav- 
ity, each with a different trap frequency determined by 
the local beam waist size. Also, Q c s can be tuned to 
as small as ~ 10 3 without significant loss in sensitivity. 
Such an instrument can thus scan over a wide variety 
of frequencies at the same time, which is crucial for the 
search for the QCD axion signals. 

Pushing the sensitivity limit in this uncharted high fre- 
quency region adds to the development of GW astronomy 
that will take place in the next decade with Advanced 
LIGO and atom interferometry (28|. GWs propagate un- 
perturbed after they are created, allowing us to study 
the remote corners of our Universe, and will become an 
indispensable tool for astrophysics and cosmology. 

We thank Savas Dimopoulos, Sergei Dubovsky, Ne- 
manja Kaloper, and Jonathan Weinstein for valuable dis- 
cussions. 
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